Halogen-Bonding-Mediated Radical Reactions: The Unexpected Behavior of Piperazine-Based Dithiooxamide Ligands in the Presence of Diiodine

N,N′-Dialkylpiperazine-2,3-dithiones (R2pipdt) were recognized as a class of hexa-atomic cyclic dithiooxamide ligands with peculiar charge-transfer donor properties toward soft electron-acceptors such as noble metal cations and diiodine. The latter interaction is nowadays better described as halogen bonding. In the reaction with diiodine, R2pipdt unexpectedly provides the corresponding triiodide salts, differently from the other dithiooxamides, which instead typically achieve ligand·nI2 halogen-bonded adducts. In this paper, we report a combined experimental and theoretical study that allows elucidation of the nature of the cited products and the reasons behind the unpredictable behavior of these ligands. Specifically, low-temperature single-crystal X-ray diffraction measurements on a series of synthetically obtained R2pipdt (R = Me, iPr, Bz)/I3 salts, complemented by neutron diffraction experiments, were able to experimentally highlight the formation of [R2pipdtH]+ cations with a −S–H bond on the dithionic moiety. Differently, with R = Ph, a benzothiazolylium cation, resulting from an intramolecular condensation reaction of the ligand, is obtained. Based on density functional theory (DFT) calculations, a reasonable reaction mechanism where diiodine plays the fundamental role of promoting a halogen-bonding-mediated radical reaction has been proposed. In addition, the comparison of combined experimental and computational results with the corresponding reactions of N,N′-dialkylperhydrodiazepine-2,3-dithione (R2dazdt, a hepta-atomic cyclic dithiooxamide), which provide neutral halogen-bonded adducts, pointed out that the difference in the torsion angle of the free ligands represents the structural key factor in determining the different reactivities of the two systems.


INTRODUCTION
Dithiooxamides are well-known organic ligands largely employed in coordination chemistry for a long time thanks to the peculiar electronic properties and coordination versatility of the two vicinal thioamidic moieties. 1−4 Despite the variety of interactions they can give rise to with metallic and nonmetallic species, dithiooxamides typically work as soft S-donors toward acceptors such as soft metal cations and halogens. 5−7 The acceptor capability of the electron density of halogens from a neutral or anionic Lewis base is well recognized as a noncovalent attractive interaction involving a halogen as an acceptor to form adducts. 8 This capability has also been elucidated by quantum mechanical methods giving rise to the widely accepted σ-hole model. 9,10 In several cases, these adducts are stable in solution and can be isolated also in the solid state, while elsewhere the system further evolves to give redox products. The driving force of further evolution of the adduct is related to the σ electronic density distribution over the three centers S−I−I and also to the stability of the formed products. 11 Numerous examples of organic substrates forming halogen-bonded adducts with iodine through their thiocarbonyl groups are known. 12 −15 Among them, N,N′dialkylperhydrodiazepine-2,3-dithione (R 2 dazdt; Chart 1), a class of hepta-atomic cyclic dithiooxamides, provides 1:1 and 1:2 diiodine halogen-bonded adducts, similarly to those obtained with other acyclic dithiooxamide donors. 8 These adducts, which have been isolated at the solid state and structurally characterized, behave as powerful lixiviants toward noble metals (Cu, Pd, Au) but are not effective toward platinum. 8,16,17 The use of these lixiviants for noble metal leaching and recovery purposes has been extensively explored for smart applications through a green chemistry ap-proach. 17 −21 On the other side, unexpected different behavior in the reaction with iodine was found for N,N′-dialkylpiperazine-2,3-dithione (R 2 pipdt; Chart 1), a corresponding class of cyclic dithiooxamides where the thioamide moieties are embedded in a hexa-atomic heterocycle. 8 Indeed, the reaction of R 2 pipdt (R = Me) with I 2 did not provide the expected halogen-bonded adduct as the final product. The X-ray diffraction (XRD) measurements at room temperature of the isolated product, dating back to the late nineties, in agreement with other analytical and spectroscopic (UV−vis−infrared (IR) and Raman) characterization, suggested the formation of the [Me 2 pipdt]I 3 , the triiodide salt of the ligand in the form of a cation. Based on the highly symmetrical configuration found for the ligand moiety by XRD (refined in the C2/c spatial group) as well as the absence of any diagnostic evidence of protonation by the other techniques, the most reasonable description of the system involved the formation of the radical [Me 2 pipdt] •+ . However, no experimental support was found to demonstrate the radical formation. This reopened the discussion on a possible protonation of the ligand. 8 Assumed this latter hypothesis, the so-called [Me 2 pipdtH]I 3 salt, differently from the previously cited adducts of R 2 dazdt, showed an exceptional reactivity toward platinum, producing the [Pt(Me 2 pipdt) 2 ]I 6 complex through a one-pot reaction in mild conditions. 22 Further studies have been performed to better characterize the [R 2 pipdtH]I 3 salts as well as to investigate the structure−property relationship of the systems based on R 2 pipdt and R 2 dazdt, in terms of both ligand interaction with halogens and leaching properties of the obtained products.
In this context, herein we report the results of the extensive studies performed over several years to provide support to the formation of the [Me 2 pipdtH] + cation and the reason why Me 2 pipdt behaves so differently from R 2 dazdt toward diiodine. Specifically, the deep characterization of a series of I 3 − salts, synthetically obtained by reacting R 2 pipdt (R = Me, Me-1; i Pr, i Pr-1; Bz, Bz-1; Ph, Ph-1) ligands with diiodine in an organic solvent, was performed. The characterization at the solid state (mainly by single-crystal X-ray and neutron diffraction measurements) and in solution (electrospray ionization mass spectrometry (ESI-MS), UV−vis), well supported by density functional theory (DFT) calculations, was able to exclude the formation of a radical species while definitely pointing out the formation of the protonated [R 2 pipdtH] + cations for Me-1, i Pr-1, and Bz-1 ligands. The finding of a different cation in the case of the Ph-1 ligand, the Ph-1C + cation resulting from an intramolecular condensation reaction, was also pointed out. A reasonable reaction mechanism for the formation of these salts, supported by DFT calculations, is also proposed and discussed to highlight why, in the case of the reaction of R 2 dazdt with diiodine, the suggested pathway does not occur.

RESULTS AND DISCUSSION
As anticipated in Section 1 and highlighted in a previous paper of ours, X-ray diffraction measurements on the product, obtained by reaction between the Me 2 pipdt ligand (Me- 1) with I 2 in organic solvents, supported the formation of a triiodide salt, [Me 2 pipdt]I 3 (Me-1I 3 ), where the cationic moiety showed a highly symmetric [Me 2 pipdt] + molecular structure. 22 However, no evidence was found to support the radical nature of this cation, which should hold an odd number of electrons ([Me 2 pipdt] •+ = Me-1 •+ ). This stimulated further studies to achieve a more reliable characterization of the salt. Moreover, open-shell DFT calculations ruled out the involvement of a radical cationic in Me-1I 3 (see Section 2.1).
The first evidence of a different cationic structure came from the ESI-MS measurements of the salt. Indeed, besides the detection of the triiodide anion (m/z 381 peak), the possible presence of a protonated [Me 2 pipdtH] + (Me-1H + ) species was supposedly based on the peak at m/z 175 ( Figure 1).  Figure S1) and the experimental details (Section S1).

Inorganic Chemistry
pubs.acs.org/IC Article 1 H NMR spectra of Me-1 (Section S4 and Figure S7) exhibit two peaks at 3.75 and 3.57 ppm corresponding to the methylene and methyl protons, respectively. Both resonances are slightly downfield-shifted in the NMR spectrum of the salt Me-1I 3 , which in addition exhibits a broad band at 2.75 ppm. Since no broadening or large shifts of the peaks, as predictable for the presence of a radical ion, were observed, a possible involvement of a radical cationic species was ruled out. While the proton signal of the S−H group is missing, the broad band at the 2.75 ppm band may be indicative of a thiol−water exchange and supported by the known typical weakness of the S−H bond. 23,24 Looking back at the original structure of the [Me 2 pipdt]I 3 salt, 22 we inferred that the existence of the Me-1H + species does not disagree with the experimental structural data. As a matter of fact, the salt crystallizes in the symmetric C2/c space group and a C2 axis is passing just in the middle of the cation, making the two parts of the molecule identical. Considering the small scattering power of the hydrogen atom, the small asymmetry in Me-1H + generated by the H + could be easily overlooked in the diffraction experiment.
An alternative way of resolving this quandary was to decrease the symmetry in the cation, substituting the methyl groups with bulkier residues, namely, isopropyl ( i Pr 2 pipdt, i Pr-   1), and phenyl (Ph 2 pipdt, Ph-1) groups. These ligands were prepared in agreement with ref 25 and, as previously made with Me-1, reacted with diiodine in a 1:2 molar ratio in CHCl 3 at room temperature (see details in Section 4). As a result, by reacting the selected R 2 pipdt with diiodine in an organic solvent, the triiodide salt of the cationic ligand was invariably obtained at the solid state and fully characterized (see Section S3, Supporting Information). Figure 2a,b shows that for i Pr-1 and Bz-1 ligands, the corresponding I 3 − salts crystallized in the P1̅ space group, hence without the C2 axis and with the asymmetric unit comprising the whole cation, differently from the previously cited Me-1. Indeed, in these two cases, a residual electron density was found away from one sulfur atom (S−H distances 1.22(3) and 1.60(6) Å and with C−S−H angles 90.5 (15) and 87(2)°, respectively, for i Pr-1 and Bz-1). Thus, a hydrogen atom was reasonably associated with such an electron density. Selected bonds and angles are reported in Table 1.
A further attempt to ascertain the presence of the thiolic hydrogen atom in the crystal structure of Me-1I 3 was made by collecting single-crystal X-ray diffraction data at low temperature (100 K) to improve the intensity and resolution of the diffraction data and reducing the thermal motion, thus allowing the identification of subtle features in the electron density. The redetermination of the crystal structure of Me-1I 3 led to a change in the space group symmetry from that previously reported C2/c to the new P2 1 /c. The discrepancy between the two symmetry group assignments was attributed to overlooked weak reflections rather than to a disorder−order transition. The new model showed no diffraction peak violations and was refined with the overall atom connectivity and molecular symmetry in agreement with the presence of ordered hydrogen on one sulfur atom. All of the hydrogen atoms, including the proton of the S−H group, were localized in the difference Fourier map calculated from X-ray data at low temperature, as shown in Figure 2c.
The nature of the bonding in the dithio-piperazinium cation Me-1H + can be discussed in detail based on its structural features in terms of bond lengths and bond angles, as described below. A similar analysis can be performed for the other presented structures that exhibit similar molecular and packing features.
The planarity of two thioamide moieties is indicative of an sp 2  Another important aspect regards the conformation adopted by the piperazine ring that is strictly connected with the dihedral angle θ between the planes comprising the thioamide moieties (N−C−S) and with the torsional angles (N−C−C− N) along the ring. The substantially smaller dihedral angles θ and the torsion angles around the central C−C bonds in the dithio-piperazinium cation (9.80, 10.89°) compared to those observed in the corresponding neutral ligand (35.28 and 34.36°) could be reasonably due to the presence of orienting intramolecular and intermolecular interactions involving the thiolic hydrogen and the triiodide ion.
The main intermolecular interactions in the crystal structure of Me-1I 3 were analyzed using the Hirshfeld surface mapped with d norm 28 and the corresponding two-dimensional fingerprint plots ( Figure 3). 29 The analysis shows that the dominant interactions in the packing are H···H, H···I, and H···S. The H··· I interactions are represented by the upper short spike (labeled a in Figure 3), in the bottom left area of the fingerprint plot. This spike is not only due to the short contact of the central iodine of I 3 − with the thiolic hydrogen atom (2.98(3) Å) but also due to multiple interactions between the lateral atoms of triiodide with the C−H groups of the ring and those of the methyl substituent on the amidic nitrogen. The lower pair of short spikes labeled b in Figure 3 represents the H···S intermolecular interactions. The spikes are due to several contacts of different lengths, involving both sulfur atoms, that overlap in the pattern.

Inorganic Chemistry pubs.acs.org/IC Article
In the fingerprint plots, many points reaching down to (1.2, 1.2) along the plot diagonal arise from the shortest H···H contacts (labeled c in Figure 3). These contacts give a strange broad diffuse region of blue points ranging from 1.2 to 1.6 Å due to multiple contacts between the C−H hydrogen atoms of adjacent molecules.
Although the positions of the thiolic protons were reliable enough to safely rule out other possible models, we decided to perform a single-crystal neutron diffraction study to obtain a more accurate determination of the hydrogen sites.
However, also in this case, the data were not good enough to allow a complete refinement of the structural parameters, as described in Section 4. Nevertheless, the crystal data obtained by these last measurements, Me-1I 3 n, were good enough to lead to a reliable position of the hydrogen and refine the crucial data: the thione moiety is confirmed to be protonated at the sulfur atom with the S−H bond distance (1.34(3) Å) in good agreement with the other values from neutron data reported in the literature. 30 A further contribution to understanding these systems came from the results of the reaction between the ligand Ph-1 and diiodine under the same conditions reported for the other compounds. Indeed, instead of the expected [Ph 2 pipdtH] + , a different structure for the cation, namely, Ph-1C + , possibly deriving from an intramolecular condensation reaction, was found, as shown in Figure 4 and detailed in Section 4.1.
As shown by the molecular structure, one sulfur atom is bonded with the ortho carbon atom of the neighboring phenyl group forming a new five-membered S,N-heterocycle in a benzothiazolylium cationic moiety. A radical mechanism mediated by iodine can be reasonably invoked for the reactions of R-1 with I 2 , as confirmed by experiments with the radical trapping (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO) (see Section S7, Supporting Information).
These results stimulated further studies on the formation mechanism of the R-1(H)I 3 salts and a comparison with the corresponding reaction undergone by R-2 ligands.

DFT Calculations and Reaction Mechanism Hypothesis.
Results from open-shell DFT calculations have contributed to ruling out the possible involvement of a radical cationic species Me-1 •+ in Me-1I 3 . The computational analysis on a potential radical Me-1 •+ species did not accurately reproduce the available X-ray experimental structure since, despite a good agreement for the S1−C1−C2−S2 dihedral angle, the calculated S1−S2 distance seems to be particularly shorter (2.87 Å) compared to the available experimental one (3.134(4) Å) in the original reported structure (Refcode QIPYIX), as shown in Figure 5. 22 The calculated S−S distance was indeed found in the range of a half-bond as it was highlighted by some of us in previous works 31−33 and close to the value of 2.8168 (11) reported for the half-bond S−S distance in the radical cation of 1,8-chalcogen naphthalenes. 34 On the other side, the experimental S−S distance larger than 3 Å can be considered nonbonding.
The occurrence of the half-bond in the calculated structure of Me-1 •+ can be explained by looking at the highest occupied molecular orbital (HOMO) (shown in Section S2, Supporting Information) of Me-1 characterized by an S1···S2 antibonding interaction. The removal of one electron from the HOMO, which generates the SOMO of Me-1 •+ (Figure 5b), weakens the antibonding character, allowing the formation of the S−S half-bond (see Figure S2 for calculated frontier orbitals). Such a behavior is also confirmed by the optimization of the Me-1 2+ species, obtained upon the removal of a second electron, generating the Me-1 2+ species, where a single S−S bond is found (calculated distance 2.22 Å; see Figure S3). On the other side, a very close agreement between the calculated (3.24 Å) and experimental (3.134 Å) S1···S2 distances, as well as between the S1−C1−C2−S2 dihedral angles, was obtained for the closed-shell [Me 2 pipdtH] + species, Me-1H + , confirming the ESI-MS observation for the presence of the protonated species in the solid state. We also excluded the protonation of one of the nitrogen atoms, the N-protonated isomer being 21.9 kcal mol −1 (91.8 kJ mol −1 ) higher in energy compared to the S-protonated one (see Figure S4).
DFT calculations have also been employed to correlate the observed different reactivity of R-1 and R-2 ligands, only differing by a −CH 2 − in the cycle, with their different structural conformation, and to elaborate a reasonable reaction mechanism. Accordingly, the reaction pathway for the general R-n with n = 1 or 2 in CHCl 3 solvent was obtained ( Figure 6). In the present discussion, only the structural data from the optimized structures are used.
In outlining a possible reaction pathway for the transformation of R-n ligands in the corresponding triiodide salts, it seems reasonable to invoke −S···I 2 halogen-bonded species of D···I 2 type. As anticipated in Section 1, depending on the donor−diiodine interaction strength, the electronic distribution changes compared to that in the free reagents, and a lengthening of the I−I distance proportional to the shortening of the intermolecular distance was observed. 11 In the presence of an excess of I 2 , another free I 2 molecule may interact with the other S atom or, as demonstrated by previous studies, may  Table 1). Thermal ellipsoids are shown at a 50% probability level. Additional data are provided in Section S5, Supporting Information. Inorganic Chemistry pubs.acs.org/IC Article promote the heterolytic I−I cleavage of the already involved I 2 molecule. The net result is the formation of halogen-bonded species SI + ···I 3 − with the SI + acting as halogen-bonding donor and the corresponding triiodide salt. 35−38 The optimized structures in the CHCl 3 solution pointed out that the energy of the two possible isomers, alternatively involving two or one sulfur center(s), differs by less than 1 kcal mol −1 (4 kJ mol −1 ) in both cases. Based on these results, it is reasonable that both isomers are present in the solution. For our purposes, we focused on the role of the second I 2 molecule as the cleavaging species for the I−I bond. In that case, the second I 2 molecule interacts through a halogen bond with the terminal I-atoms of the chain. The following step is the heterolytic cleavage of the I1−I2 bond, resulting in an I 3 − anion and a ligand−I + cation, as well demonstrated in the case of 1,3-bis(2,6-diisopropylphenyl)-1,3-dihydro-2H-imidazole-2thione in the presence of one and two iodine molecules. 39 Taking the above as typical behavior of soft donor−acceptor systems, it has been applied to both the cited classes of ligands, specifically indicated as Me-1 and Me-2.
As shown in Figure 6, starting from the free ligand, the reaction pathway consists of five main steps. It is worth noting by comparing the two ligands that the main structural difference for Me-1 and Me-2 is in the value of the dihedral S−C−C−S angle calculated as 40 and 74°, respectively. In the first step, the formation of an L−I 2 halogen-bonded adduct is considered. The main geometrical difference between Me-1·I 2 and Me-2·I 2 remains in the S−C−C−S angle (45 vs 70°for Me-1·I 2 and Me-2·I 2 , respectively). Furthermore, the I1−I2 bonds are partially elongated compared to the free I 2 molecule (e.g., 2.95 Å Me-1·I 2 against 2.74 Å), while the S1−I1 distance is around 2.87 and 2.85 Å for Me-1·I 2 and Me-2·I 2 , respectively. The next step is the interaction of a second I 2 unit with the coordinated one, which causes a weakening of the I1−I2 bond and the strengthening of the S1−I1 bond. As a Figure 6. Proposed reaction mechanism for the formation of Me-1H + and Ph-1C + species and comparison with the corresponding Me-2 system based on the calculated structures in the CHCl 3 solvent.
Step 3 reports the calculated structures of the Me-1I + and Me-2I + , where the I 3 − group leaves the molecular compound. For Me-2I + , the S1−I1 distance is 2.48 Å and the C1−S1−I1 angle is 107°(with an S−C−C−S dihedral angle of 71°). Calculated values for I1−I2 distances in Me-1·(I 2 ) 2 and Me-2·(I 2 ) 2 and S1−I1 in Me-2I + well agree with literature crystal data for other similar species (see the corresponding specie obtained by 1,3-bis(2,6-diisopropylphenyl)-1,3-dihydro-2H-imidazole-2-thione, as a reference). 39 On the contrary, Me-1I + presents several unusual features. First, the S1−I1 distance is 2.60 Å, 0.12 Å longer than in Me-2I + and in a typical S−I single bond. More interestingly, the S1−S2 distance is as short as 2.78 Å, a distance typical of a half S−S bond. For comparison, the S1−S2 distance in Me-2I + is longer than 3.40 Å. Moreover, the S1−C1−C2−S2 dihedral angle is almost flat (7.3°) and the S2−S1−I1 angle is 155°, close to the linearity. The geometrical arrangement in Me-1I + is in agreement with a three-center-four-electron hypervalent system description. 40 Accordingly, in the case of Me-1, the presence of the hypervalent bond suggests the irreversible homolytic splitting of Me-1I + in two radicals, namely, Me-1 •+ and I • . Remarkably, as highlighted in Figure 8, the energy required for this splitting is less than half for Me-1I + with respect to Me-2I + .
The calculated free-energy profiles, shown in Figure 9, point out the significantly more favored formation of Me-1I + than the corresponding Me-2I + species, attributable to the presence of the hypervalent bond in Me-1I + .
In the final step, the recombination of the radical species is considered. Specifically, Me-1 •+ will compete with an R−H for an H • extraction to provide the Me-1H + structure. The most plausible R−H source of the radical H • is the CHCl 3 solvent. The same behavior is expected for similar systems based on i Pr-and Bz-substituents. In the case of the phenyl-substituted Ph-1, Ph-1 •+ will go for an internal condensation reaction to produce the radical intermediate Ph-1CH •+ . The latter, featuring a tetrahedral carbon center, may interact with radical species in solution to provide the final species Ph-1C + together with some radical recombination products. It is worth noting that the Ph-1CI 3 condensation product is an indirect proof of the existence of an R-1 •+ specie in solution where the phenyl group is acting as an internal trapping radical. Besides, I • may combine with R • and/or with I − or I 3 − to form I 2 •− , followed by a disproportion reaction to I 3 − and I − , in agreement with the previous finding by other authors. 43, 44 To summarize, despite the possibility for the two systems to evolve toward the R-N·2I 2 (N = 1, 2) adducts, in the case of R-1 ligands, the irreversible formation of an R-1 •+ radical is   Inorganic Chemistry pubs.acs.org/IC Article energetically feasible shifting the reaction toward the formation of the triiodide salts. On the contrary, for Me-2 only an equilibrium between the free ligand and various iodine adducts is envisaged due to the energetically disfavored formation of the radical species. As a consequence of the above, we can state that the difference in the torsion angle of the free ligands has been proved to represent the structural key factor in determining the reactivity of the two systems.

CONCLUSIONS
The reaction of the hexa-atomic cyclic dithiooxamide ligands (R 2 pipdt, R-1) with diiodine, unexpectedly providing triiodide salts, has been thoroughly revisited and newly investigated. Extensive experimental work and new theoretical interpretative frameworks were required for shedding light on the structural features of the reaction products and the reaction mechanism behind the triiodide salt formation. At the end of this journey, X-ray diffraction measurements, complemented by neutron diffraction experiments, were able to experimentally establish the main formation of S-protonated ligands upon the reaction of R-1 ligands (with R = alkyl and Bz) with iodine in organic solvents. Differently, a new five-membered S,N-heterocycle in a benzothiazolylium cationic moiety was formed in the case of R = Ph. DFT calculations, based on the structural data of Me-1 and Me-2, allowed us to speculate on a reasonable reaction mechanism. Specifically, when halogen bonding occurs with I 2 in the case of Me-1, the almost planarity of the S−C−C−S moiety allows the formation of a three-center-four-electron S− S−I bond, which irreversibly evolves toward active radical species through the homolytic cleavage of the S−I bond. Eventually, the radical ligand R-1 • induces an intra-(as for R = Ph) or inter-(as for the other R substituents) molecular radical reaction with the formation of the corresponding I 3 − salts. Otherwise, in the case of the Me-2, the wider S−C−C−S torsion angle does not allow the formation of the three-centerfour-electron S−S−I bond hindering the radical activation with the consequence that only Me-2·nI 2 halogen-bonded adducts occur.
Thus, the reported combined experimental and computational study allowed us to elaborate on a reasonable reaction mechanism providing an interesting example of a halogenbonding radical reaction promoted by diiodine. This opens the way to further investigations on the peculiar reactivity demonstrated by R-1HI 3 salts toward metals, primarily platinum.

Materials and Methods.
Reagents and solvents of reagentgrade quality were used as supplied by Sigma-Aldrich. The R 2 pipdt ligands R-1 were prepared according to ref 25. Ph 2 pipdt was synthesized as described in the following section.
SIR2004 48 and SHELXL programs 49 were used for structure solution and refinement on F 2 by full-matrix least-squares techniques. All non-hydrogen atoms were refined using anisotropic displacement parameters. The H atoms were located from different Fourier syntheses ( i Pr-1I 3 and Me-1I 3 ) or placed in geometrically calculated  Inorganic Chemistry pubs.acs.org/IC Article positions (Bz-1I 3 , Ph-1CI 3 , and Ph-1) and included in the refinement using a riding model in conjunction with a U iso (H) = 1.5 U eq (CH 3 , SH) or U iso (H) = 1.2 U eq (CH 2 ) constraint. The diagram was drawn using ORTEPIII program. 50 Crystal data and structure determination results are summarized in Table 2.

Neutron Diffraction.
Neutron diffraction experiments on Me-1I 3 complexes (Me-1I 3 n) were carried out using the very-intense vertical-axis Laue diffractometer (VIVALDI) 51,52 at the Institut Laue-Langevin, Grenoble, with a white neutron beam covering wavelengths from 0.8 to 5.2 Å. The crystal (dimensions) was glued on a vanadium pin using an epoxy resin, mounted on the diffractometer, and cooled down to 100 K in an orange He flow cryostat.
To ensure full coverage of reciprocal space, a total of 11 Laue diffraction patterns were collected at 15°intervals in rotation about the vertical axis perpendicular to the incident beam, with an exposure time for each frame of around 3.8 h.
The Laue patterns were indexed and integrated using the LAUEGEN software 53 suite, and wavelength normalization was carried out using the LAUENORM program. 54 Correction for absorption was deemed unnecessary due to the small crystal dimensions.
Since the indexing of a Laue pattern provides only relative unit cell dimension, the absolute unit cell parameters were determined from Xray, working at the same temperature as the neutron analysis. Data were refined by full-matrix least-squares starting from existing X-ray models based only on the atomic coordinates for the heavy atoms while all hydrogen atoms were located from Fourier difference maps.
Considering very limited sample dimensions, these collected neutron data were of sufficiently good quality to find hydrogen positions and confirm the S−H assignment. However, only individual isotropic factors were used in the refinement.
The collected data were of poor quality but good enough to find hydrogen positions and confirm the S−H assignment. However, only individual isotropic factors were used in the refinement.

Computational Details.
All of the compounds were optimized with the Gaussian16 suite of program 55 using the hybrid density functional B97D. 56 All of the free energies, derived after the calculations of the vibrational frequencies, refer to a temperature of 298 K. All of the calculations were based on the CPCM model 57, 58 for the chloroform solvent. The basis set 6-31G inclusive of polarization functions was used for all species, while for iodine the Stuttgart/ Dresden (SDD) pseudo-potential 59 was employed. The coordinates of the optimized structures and their energetic parameters are reported in Section S6, Supporting Information. ■ ASSOCIATED CONTENT